Plasmons do not go that quantum 
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We develop a theoretical model of the surface plasmon resonance of metallic nanospheres in the 
size range down to the single nanometer size. Within this model we explicitly show how different 
microscopic mechanisms, namely quantization due to size (QSE) and electron spill-out, affect the 
energy of the surface plasmon. We demonstrate, that electron spill-out effects, which can move 
the surface plasma energy both toward the red or the blue, can be comparable to or even stronger 
than QSE. Thus, depending on circumstances, QSE may only be observed for ultrasmall metal 
nanoparticles much closer to 1 nm in size than to 10 nm. Results presented herein are in quantitative 
agreement with recent published experimental results for Ag and Au. 
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During the last 20 years there has been considerable 
progress in both calculating and measuring the properties 
of small particles 1 - - — as well as the corresponding efforts 
for clusters. Of particular interest is to monitor the re- 
gion where cluster research meets and merges with the 
corresponding small particle research^—. Furthermore, 
in the last few years, the behavior of plasmon resonances 
in nanosystems has become a hot topic for its important 
applications in varied fields ranging from cancer therapy, 
nanophotonic devices, biosensing to catalysis. Interest- 
ing physical aspects of these resonances include, among 
others, investigation of the quantum regime of tunnel- 
ing plasmonics^ or what is an ultimate limit to electro- 
magnetic enhancement in these systems^ 1 -" . In spite of 
all this effort, several aspects of the behavior of plas- 
mon resonances in nanometer-sized metallic spheres are 
not fully understood. A particularly interesting exam- 
ple concerns the energy shift experienced by the surface 
plasmon as a function of the size and shape of the parti- 
cle. There has been in the literature a variety of experi- 
mental reports of blueshifts^^r— or redshifts^-^ of these 
resonances as the particle dimensions arc reduced, which 
today may seem conflicting. There have, of course, been 
many attempts during the years to address this question 
theoreticall y 14 ' 20 7 ~— . However, it is only recently that 
the experimental situation has become so well controlled 
that a realistic comparison can be made with theoretical 
model predictions. 

The classical Mie electromagnetic theory 2 ^ is adequate 
for describing the optical properties of particles with a 
radius larger than 10 nm. Its main ingredient is the fre- 
quency dependent dielectric function e(ui) of the material 
and its most characteristic feature is the presence of a sur- 
face plasmon resonance. The resonance moves in energy 
as the particle size is changed, depending on the value 
of the radius R with respect to the effective wavelength 
of the plasmon c/w p , where c is the speed of light and 
Lo p is the plasma frequency of the volume plasmon. The 
Mie theory gives rise to saturation of this energy at the 
value Lof = uip/s/1 + 2e m for radii smaller than the wave- 



length, where e m is the dielectric constant of the medium 
surrounding the particle (assumed to be frequency inde- 
pendent). When dealing with particles in the size-range 
1-10 nm the surface starts to become mor e and more 
important compared with the bulk response and thus in- 
fluences heavily collective modes localized in the surface 
region. Since electrons are able to spill out of the metal 
because of the finite potential barrier at the surface, the 
surface screening of an external perturbation is drasti- 
cally changed from the Mie model assuming, as it does, 
a sharp interface between a particle and the surrounding 
medium. Whereas the classical theory is characterized 
by bulk dielectric properties and whence the length scale 
R (the particle radius), the surface enters the electrody- 
namic response as a frequency dependent complex length 
scale d r (w)2£. Its real part is a measure of the center of 
gravity of the screening charge induced in the particle by 
an external perturbation of frequency to and its imagi- 
nary part describes surface absorption. Reducing the ra- 
dius below a few nanometers we move to a region where 
again the bulk becomes important i n a new fundamen- 
tal way entering the quantum size effect region where 
the discrete energy level spectrum makes it possible to 
set up standing electron-hole pairs and quantum core 
plasmons^. The last are charge oscillations localized 
at the very center of the particle because of the quan- 
tum confinement of the electrons. These quantum size 
effects influence the particle response to a significant de- 
gree and compete with the effect coming from the surface. 
For these ultra small particles sophisticated Time De- 
pendent Density Functional Theory (TDDFT) and Time 
Dependent Local Density Approximation (TDLDA) cal- 
culations are available nowadays 1 ^^^. However, these 
become increasingly computationally demanding with in- 
creasing radius. 

The purpose of this article is to clarify the role of the 
different microscopic mechanisms affecting the energy of 
the surface plasmon of metallic spheres as a function of 
their size, for diameters going down to a few nanome- 
ters. To do so we develop a theoretical model including 



surface spill-out and quantum size effects. This model is 
based on known concepts of surface screening and quan- 
tum size effects and successfully bridges cluster (sub 1 
nanometer) research and Mie theory (greater than 10 
nm) predicting results that agree with TDDFT calcu- 
lations. Our work is triggered by recently available very 
detailed cxpcrimcnt o 17 i 18 which greatly improve the pos- 
sibilities of checking basic electromagnetic response the- 
ory predictions with experimental findings. We note that 
the theoretical framework used in these papers does not 
appropriately include the major effect of the surface spill- 
out electrons thus giving the reader the impressi on that 
quantum effects start to enter the picture for diameters 
much larger than they actually do. We find that electron 
spill-out effects, which can move the surface plasma en- 
ergy both toward the red or the blue, can be comparable 
to or even stronger than QSE. We unambiguously show 
that taking surface screening and d-electrons properly 
into account the experimental surface plasmon energies 
of Ag and Au particles in a variety of host media can 
be reproduced down to a few nm sized particles and the 
discrepancy at smaller sizes can be accounted for in a 
simple model of quantum size effects. 

The electromagnetic response of a particle in the dipo- 
lar approximation can be described by its polarizability. 
The polarizability of a metal sphere of radius R with a 
bulk dielectric constant e(u>) can be written as2£r— 



a(w) = R s 



( £ H- Cm )(l-*M) 
e(w) + 2e m + 2 (e(w) - e m ) 



d T {u) 
R 



(1) 



where the length d r (u>) is related to the charge density, 
5p(r,Lo), induced by the external perturbation, by the 
formula 



d r {u>) J dr r(R — r)5p(r,uj) 
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where r is the radial coordinate and the integrals ex- 
tend to the whole space. Notice, that for the Mie model, 
where the induced charge density is right at r — R, d r is 
identically zero and eq. ([1} reduces to the classical Mie 
result when retardation is neglected. We will focus in 
this paper on the shift in energy of the surface plasma 
resonance, which can be obtained from the poles of the 
polarizability. Then, from eq. ([1]), the frequency of the 
surface plasmon resonance, u) s , has to fulfill 
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It has been known long since22r— that the calculated 
energies of the surface plasmons are very dependent on 
the shape of the barrier that confines the electronic sys- 
tem, even for semi-infinite metals. Infinite barrier models 
confining the electrons inside the metal produce frequen- 
cies higher than the classical value while realistic surface 
barriers can produce lower frequencies because electrons 



are allowed to spill-out of the metal, thus lowering the 
plasma frequency. The same happens for nanospheres. 
If electrons are allowed to spill-out of the metal, the 
induced surface charge density has in general its main 
weight outside the surface with its center of gravity giv- 
ing Re[d r (w)] < 0, while a positive value of Re[d r (w)] is 
obtained if electrons spill-in. From eq. ([3]) the color of 
the shift of the surface plasmon frequency with respect 
to oof is directly related to the sign of Re[d r (w)]: blue 
for positive and red for negative signs, respectively. Con- 
sequently, theoretical models in which the electrons in 
the nanoparticle are confined directly or indirectly by an 
infinite barrie r 17 ' 18 i 32 can only produce blueshifts. 

The above considerations refer to surface effects. Now 
we analyze quantum size effects (QSE) which are also 
present at the nanoscale size of interest here. These ef- 
fects can be described approximately by a dielectric con- 
stant similar to the one of a semiconductor: the discrete 
nature of the electronic states requires a minimum of en- 
ergy to excite an electron, which is equivalent to having 
an energy gap. Inspired by the Penn dielectric function 
for semiconductors 3 ^, we propose to use the following 
form of the real part of the dielectric function: 



e(w) = 1 - 



(4) 



where, as introduced by Gorkov and Eliashberg^i, the 
energy gap A scales with the particle radius as A = 
uipRo/R. Then, using a simple box model, one obtains 
from their treatment Rq ~ Llaoy^, where r s is the ef- 
fective free-electron density parameter and ao is the Bohr 
radius. The existence of an energy gap influences the col- 
lective modes of the sphere. The classical bulk plasma 
frequency given by e(w) = 0, depends on the radius R as 



v P (R) 




(5) 



and, in the absence of any surface effect, the frequency 
of the surface plasmon changes to 



u s (R) 




(6) 



That is, QSE will produce a blucshift of both surface and 
bulk plasmon energies. We notice that the relative shift 
is larger for the surface compared to the bulk plasmon. 

Consequently, in a real situation with both quantum 
and surface barrier effects present, there is a competition 
between the blue- and the redshifts. Hence, in this work 
we apply the two just outlined models to calculate surface 
plasmon energies in small particles which will be com- 
pared to experimental results published in Refs i 16 i 17 ' 35 . 
Since we will be dealing with the noble metals Ag and 
Au, we also need to include the screening due to the d- 
electrons by modifying eq. dU as 



e(oj) = e d (u>) 
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FIG. 1: Surface plasmon energy versus diameter for small Ag 
nanospheres. Points with error bars are experimentally ex- 
tracted results from EELS measurements by Scholl et al.— . 
The blue line is our calculation of the surface plasmon energy 
shift resulting from electron spill-out and the magenta line 
includes also the effects of size quantization. It is clear that 
the blue line gives a reasonable account for the data but the 
magenta line shows an even better agreement with the exper- 
iment. It has earlier been claimed that quantum effects start 
already at 10 nm however our calculations show that they do 
not dominate until the particles are about 10 times smaller. 



where ed(w) is obtained from experimental optical data£& 
This form of e(w) is then introduced into eq. ^ to take 
into account quantum size effects and d-electrons on the 
surface plasmon frequency and we find that it has to fulfill 
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In the detailed treatment above we looked at the center 
of gravity of the induced density d r . We now ask what 
is the influence of d-electrons and quantum size effects 
directly on d r . We can do this to lowest order by looking 
to the corresponding quantity for a planar surface d± (w) 
because it has been shown^ that the induced charge den- 
sity at the surface of a sphere is very similar to that of 
a planar surface down to a few nanometers in size. Also, 
we note that the d-electrons are very localized and there- 
fore largely excluded from the surface region where the 
conduction electrons spill-ou t 37 ' 38 . Hence, we will use 
the values of d±(uj) calculated by Feibelman 3 ^ for a pla- 
nar jellium surface of r s =3 (r s = 3.02 and r s = 3.01 
are the values giving the density of conduction electrons 
in Ag and Au, respectively). However, the screening of 
the d-e lectrons is very important in the bulk, in partic- 
ular it affects the bulk plasma frequency by changing uj p 
from its free-electron value of 9.1eV to w* = 3.81 eV and 
oj* = 5.99 eV for Ag and Au, respectively. Since d± only 
depends on the ratio lo/lu p , we use for d r the values of 
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FIG. 2: Surface plasmon energy shift for Ag particles embed- 
ded in a solid argon matrix. Full dots are the experimental 
data of Charle et al>2£. The blue line shows our results from 
electron spill-out and the magenta line includes also the ef- 
fects of size quantization. 



d± with u> rcnormalizcd to uj* . Then wc introduce the ef- 
fects of d-electrons and of size quantization on d r trough 

u;(R). 

We apply eq. $8§ to calculate the energies of the surface 
plasmon resonances of three different metals: Ag, Au, 
and a pure free-electron-metal of r s = 3 . These metals 
have in common their close values of the effective free- 
electron density parameter r s . However, the screening 
caused by the d-electrons is very different in Ag, Au, 
and a free-electron metal at the typical frequencies of 
the surface plasmons. This is a crucial physical property 
for making the resonances of metal nanoparticles to shift 
to different colors. 

Figurc[T]shows the calculated frequencies of the surface 
plasmons of Ag nanospheres as a function of the diame- 
ter, together with the experimental results by Scholl et 
al. for Ag spheres deposited on carbon foils obtained us- 
ing electron energy loss spectroscopy (EELS)^ 7 - (points 
with error bars). In our calculation we use e m = 1.5, 
as an average of the dielectric constant of the medium 
surrounding the nanospheres, which reproduces the mea- 
sured surface plasmon energies at the largest sizes. Even 
though there is large scatter in the experimental data 
there is a clear trend of increasing plasma frequency as 
the particles get smaller. The blue line is our calcula- 
tion of the surface plasmon energy shift resulting from 
electron spill-out. Usually this spill-out gives a decrease 
in surface plasmon energy but for Ag the screening of 
the d-electrons can in some circumstances push the sur- 
face plasmon energy into a region where the dynam ical 
Re[d r (w s )] is positive and gives a blueshifti^. It is clear 
that the blue line gives a reasonable account for the data 
except at the very smallest particles of a size 1-2 nm. 
In this region we do expect effects from the finite size 
of the particles, introducing discrete energy levels and a 
gap at the Fermi level. Accounting for this as indicated 
by eq. ([7]) and using the simple box estimation of Rq, 



2.8 



> 

&2.6 

i_ 
Q) 

§2.4 



2.2 




■ Ag, Lerme et al., Cottancin et al. 
O- Ag, QSE + spill-out 
■A- Ag, spill-out 

• Au, Lerme et al., Cottancin et al. 
O- Au, QSE + spill-out 
V Au, spill-out 



5 7 

diameter [nm] 



FIG. 3: Surface plasmon energy positions of Ag and Au parti- 
cles embedded in porous alumina. Full black symbols without 
a line represent the experimental data for Ag (squares) and 
Au (dots) from works of Lerme et al.— and Cottancin et al.— . 
The red lines represent our results from electron spill-out and 
the blue lines also include QSE. We see here that for Ag (top) 
we are in a redshifted spill-out region. However for very small 
particles the quantum size effects contributes a shift in the op- 
posite way which even upsets the redshift. In the bottom part 
of the figure we illustrate that our treatment also works for 
other metal particles, in this case gold. 



Rq = l.laoy 7 ^, eq. © results in the magenta line yield- 
ing an even better agreement with the measured data. 
We would like to stress that the only free parameter of 
this calculation is e m , which, however, is constrained in 
its value, since it has to yield the asymptotically mea- 
sured plasmon energies. Furthermore, the resulting value 
is very close to the average of air and carbon used in the 
experiments as support. 

The good agreement between theory and experiment 
found for Ag particles on carbon foils is not fortuitous. 
We show this by analyzing other experiments in which Ag 
and also Au nanospheres are embedded in solid matrixes. 
Figure [2] shows the calculated frequencies of the surface 
plasmons of Ag together with the experimental results 
by Charle et ali^£ for Ag spheres embedded in a solid 
argon matrix. We have used e m = 1.85 which is very 
close to the value e m = 1.75 quoted in^ for solid argon. 
The blueshift of the resonance is explained exactly on 
the same basis as the results of Fig. Q] and again this 
simple model is able to reproduce the experimental data 
on a quantitative level. Notice that the maximum value 
of the shift measured (and calculated) is approximately 
0.3 eV while it was 0.5 eV in Fig. Q] This is due to the 
fact that as e m increases, tof decreases thus moving to 
the a region of frequencies where the dynamical Re[<i r (w)] 
is still positive but smaller than in the previous case. 

Figure[3]shows the calculated frequencies of the surface 
plasmons of Ag and Au nanospheres together with the 
experimental results by Cottancin et al.— and Lerme et 
alJ£ for Ag and Au spheres embedded in alumina (black 
symbols). In all the calculations we have used the value 
e m = 2.8 appropriate for porous alumina^. As in Figs. 
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FIG. 4: Surface plasmon energy positions of free-electron- 
metal particles of r a — 3 in vacuum. The redshift produced 
by surface spill-out cannot be overcome by effects of size quan- 
tization as in the case of Au and Ag and the surface plasmon 
moves to lower energies with decreasing size. 



Q] and [2j we show results considering only surface effects 
and including effects of size quantization as well. 

First, we discuss the case of Ag. Notice that the rele- 
vant surface plasmon energies now are pushed down be- 
low 3 eV as compared to Figs. Q] and [2] due to the differ- 
ent dielectric properties of alumina as compared to car- 
bon/vacuum and solid argon, respectively. Here we find 
that surface effects cause a redshift of the resonance. This 
is because, as we have anticipated, an increased value of 
e m has moved the classical Mie frequency to the region 
of negative values of Re[d r ] thus changing the sign of 
the shift. Then for Ag we are usually in a region where 
the spill-out effect can go both red and blue for a small 
change of frequency. In the present case the effects of size 
quantization compensate the redshift in the way of keep- 
ing the energy of the surface plasmon almost constant as 
the size decreases, as seen in the experiment. 

In the case of Au, we also obtain that surface effects 
would shift the frequency toward the red but now the 
blueshift resulting from size quantization is stronger, pro- 
ducing the net blueshift shown in the figure. For Au we 
use an Rq which is 40% larger than for Ag. A possible 
reason for this is that, even though electrons in the d- 
bands will in principle also experience size effects, these 
should be smaller for Ag with deeper d-bands than for 
Au. Notice in eq. ([5]) how the relative surface versus 
quantum size effects depend on €d{^>)- This effect is il- 
lustrated in the calculations we present next, for which 
6d(w) = 1. 

Figure 3] shows our results for frce-electron-metal 
nanospheres of r s =3 in vacuum, using as before the 
simple box estimate for Rq. This is an example where the 
surface effects are practically responsible for the strong 
redshift of the surface resonance, size quantization ef- 
fects being unable to compensate for it. We would like 
to mention that sophisticated calculations for this sys- 
tem have produced redshifts of the surface resonances^. 
The value we find for a radius of R = 0.74 nm is 



uj s = 0.94:Cjf which is close to the value ui s = 0.91ujf 
obtained by Townsend and Bryant^ using TDDFT. More- 
over, a net redshift of these resonances has been measured 
for nanospheres of the typical free-electron metals Na and 
K in Ref£ and the same behavior is expected for other 
such metals as Al and Mg. 

We have developed a model theory that is in good 
agreement with existing sophisticated calculations only 
feasible for very small systems and yields excellent quan- 
titative agreement with a variety of available experimen- 
tal results for metallic spheres of nanometric size. We 
hnd that the effects of the spill-out of the electronic 
charge at the surface of the nanoparticle, which can 
move the surface plasma energy toward the red or the 
blue depending on the circumstances, are never negligi- 
ble in comparison to effects of size quantization which 



are also present in these small systems. We thus con- 
clude, that the surface plasmon resonance in small par- 
ticles does not go that quantum as recent experiments 
seem to sugges t ' 18 : instead of 10 nm we would expect 
to see effects unambiguously rather at 1 nm and below. 
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